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ABSTRACT 

We study the four BL Lac objects (RGB J0152+017, 1ES 0229+200, 1ES 0347-121 and PKS 
0548-322) detected in the TeV band but not present in the 1FGL catalogue of the Fermi/Large 
Area Telescope. We analize the 24 months of LAT data deriving 7-ray fluxes or upper limits 
that we use to assemble their spectral energy distributions (SED). We model the SEDs with a 
standard one-zone leptonic model, also including the contribution of the reprocessed radiation 
in the multi GeV band, emitted by the pairs produced through the conversion of the primary 
TeV emission by interaction with the cosmic optical-IR background. For simplicity, in the 
calculation of this component we adopt an analytical approach including some simplifying 
assumptions, in particular i) the blazar high energy emission is considered on average sta- 
ble over times of the order of 10 7 years and ii) the observer is exactly on-axis. We compare 
the physical parameters derived by the emission model with those of other high-energy emit- 
ting BL Lacs, confirming that TeV BL Lacs with a rather small GeV flux are characterized 
by extremely low values of the magnetic field and large values of the electron energies. The 
comparison between the flux in the GeV band and that expected from the reprocessed TeV 
emission allows us to confirm and strengthen the lower limit of B > 10 -15 G for the inter- 
galactic magnetic field using a theoretically motivated spectrum for the primary high-energy 
photons. 

Key words: radiation mechanisms: non-thermal — 7-rays: theory — 7-rays: observations — 
galaxies: individual: RGB J0152+017 — galaxies: individual: 1ES 0229+200 — galaxies: in- 
dividual: 1ES 0347-121— galaxies: individual: PKS 0548-322 — galaxies: individual: PMN 
J0148+0129 
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1 INTRODUCTION 

Blazars, characterized by a relativistic jet pointing towards the 
Earth, represent the most extreme flavor of Active Galactic Nuclei 
(AGN). Their relativistically-boosted non-thermal continuum, ex- 
tending from radio to 7-ray energies, is characterized by a spectral 
energy distribution (SED) with two broad peaks, with maxima at 
IR/UV or even X-ray band (due to synchrotron emission by rela- 
tivistic electrons in the jet) and in the 7-ray band (probably pro- 
duced through inverse Compton scattering), respectively. Fossati et 
al. (1998) proposed the existence of an inverse relation between the 
frequencies of both peaks and the (radio or bolometric) luminos- 
ity, the so called "blazar sequence". The most powerful sources are 
thus characterized by peaks located at low frequencies, while the 
less powerful ones display the highest peak frequencies (hence the 
classification as Highly peaked BL Lac objects, HBLs, Padovani 
& Giommi 1995) and their high-energy bump can reach TeV en- 
ergies. Sources detected at Very High Energies (VHE, E > 100 
GeV) still form a small populatior[j by they are intensively stud- 
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ied (e.g. De Angelis et al. 2008), since they represent ideal natural 
laboratories to investigate particle acceleration and cooling and to 
indirectly probe the extragalactic background light (EBL). 

Tavecchio et al. (2010a, hereafter TlOa), studying all the BL 
Lac of the LBAS sample (Abdo et al. 2009; it includes all the AGN 
detected at high significance by the Large Area Telescope onboard 
Fermi in the first three months of operation), and/or at VHE by 
Cherenkov telescopes, discussed the possible existence of two sep- 
arate groups of high-energy detected BL Lacs. In fact, while the 
majority of them, when fitted with a standard one-zone synchrotron 
self Compton model, are characterized by values of the magnetic 
field in the range between 0.1-1 G and Lorentz factors of the elec- 
trons emitting at the peak of the order of 7b ~ 10 4 — 10 5 , a hand- 
ful of sources are characterized by very low magnetic fields (of the 
order of 10~ 2 - 10~ 3 G) and very large electron energies, corre- 
sponding to Lorentz factors of 7b ~ 10 6 — 10 7 . These peculiar 
values of the physical quantities are directly connected to the ex- 
tremely high frequencies of the peaks in the SED (hence the name 
of "extreme" HBLs introduced by Costamante et al. 2001). 

Because of the extremely low flux characterizing these sources 
in the GeV band, some of these extreme HBLs have been recently 
considered for deriving upper limits on the still unknown inten- 
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sity of the intergalactic magnetic field (IGMF, e.g. Kronberg 2001, 
Widrow 2002), based on the measured upper limits of the flux in the 
GeV band (Neronov & Vovk 2010, Tavecchio et al. 2010b, here- 
after TlOb, Dolag et al. 201 1). 

Constraints on the magnetic field intensity in voids could help 
to understand the origin of the "seed" fields assumed in dynamo 
amplification models for magnetic fields in galaxies and galaxy 
clusters (e.g. Kulsrud & Zweibel 2008). Proposed models for these 
seed fields range from those considering effects occurring dur- 
ing the inflation (e.g. Turner & Widrow 1988) or phase transition 
era in the Early Universe (e.g. Kahniashvili et al. 2011 and refer- 
ences therein) to those invoking mechanisms active during the early 
stages of protogalaxy formation (e.g. Gnedin et al. 2000). Classical 
methods used to measure or derive upper limits on IGMF consider 
Faraday rotation of polarization angle of radio emission of quasars 
(e.g. Kronberg 2001) or the effects of magnetic fields on the Cos- 
mic Microwave Background (e.g. Durrer et al. 2000). 

A rather promising method to obtain lower limits to the field 
(hence complementary to the methods discussed above) is based on 
the fact that VHE photons emitted by blazars are absorbed and con- 
verted into electron-positron pairs through the interaction with the 
optical-IR cosmic background (e.g. Nikishov 1962). These pairs, in 
turn, rapidly cool through inverse Compton scattering off photons 
of the CMB, emitting photons of lower energies. If the primary 
spectrum extends up to very high energies, part of the reprocessed 
emission will still be above the threshold for absorption and a fur- 
ther generation of pairs will be created. If the maximum energy of 
the primary emission is sufficiently large this process will lead to 
the formation of an electromagnetic cascade (e.g. Aharonian et al. 
1994, Coppi & Aharonian 1997). While the pairs emit, they are de- 
viated by the original trajectory by the tiny IGMF. Therefore, the 
observed reprocessed emission will be spread over a solid angle 
larger than the original one, "diluting" its flux: the comparison be- 
tween the predicted and the observed reprocessed flux provides a 
direct measure or a lower limit of the IGMF intensity (Plaga 1995, 
Dai et al. 1992, D' Avezac et al. 2007, Murase et al. 2008, Dolag et 
al. 2009, Elyv et al. 2009). 

In this paper we consider the extreme HBLs not detected 
by LAT in the first 11 months of observations (1FGL catalogue). 
In this catalogue four TeV HBLs are not present, namely RGB 
J0152+017 (z = 0.080), 1ES 0229+200 (z = 0.14), 1ES 0347- 
121 (z = 0.188) and PKS 0548-322 (z = 0.069). 

The aim of this paper is twofold: 1) we intend to model the 
SED of the sources and derive the physical parameters consider- 
ing also the possible role of the reprocessed emission in shaping 
the observed high-energy spectrum; 2) we use the modeling of the 
primary blazar high-energy emission to calculate self-consistently 
the reprocessed component and therefore derive the lower limits on 
the IGMF. We stress that both issues are necessarily interlaced: the 
reprocessed component can provide a non-negligible contribution 
to the 10-100 GeV emission (or even above), while the primary 
TeV spectrum is necessary to have a reliable characterization of the 
expected reprocessed component. 

For all the sources we analize the ~24 months LAT data deriv- 
ing upper limits (or fluxes, in few cases) in the 0.1-100 GeV band 
that we use, together with data at other frequencies, to construct 
their SEDs (§2). We model their SEDs with the standard one-zone 
leptonic model and we derive the physical parameters of the emit- 
ting region (§3) also considering the possible contribution of the IC 
emission from the pairs produced by the conversion of high-energy 
photons interacting with the EBL, whose level and spectral shape 



Energy band [GeV] 
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TS 


0.1-1 


111 


17 


1-10 


45 


56 


10-100 


5 


43 



Table 1. Results of the analysis of the LAT data for RGB J0152+017. For 
each energy band we report N, the number of photons predicted by the 
model, and the corresponding value of the test statistics, TS 

also depends on the value of the IGMF. In §4 we show the results 
and in §5 we conclude. 

We use Ho = 70 km s _1 Mpc" 1 , 0\ = 0.7, Q.m = 0.3. 



2 SPECTRAL ENERGY DISTRIBUTIONS 
2.1 LAT data 

RGB J0152+017, 1ES 0229+200, 1ES 0347-121 and PKS 0548- 
322 do not appear in the 1FGL catalogue (Abdo et al. 2010), which 
includes all the sources detected above 100 MeV with a signicance 
exceeding 4.5a during the first 1 1 months survey of LAT onboard 
Fermi (Atwood et al. 2009). 

We used the the publicly available dat^l to search for detec- 
tions or derive updated upper limits. We selected the photons of 
class 3 (DIFFUSE) with energy in the range 0.1-100 GeV col- 
lected from 2008 August 4 (MJD 54682) to 2010 July 15 (MJD 
55392), for a total of about 24 months of elapsed time. These data 
were processed by using Science Tools 9.15.2, which includes the 
Galactic diffuse and isotropic background and the Instrument Re- 
sponse Function IRF P6 V3 DIFFUSE. 

We selected photons in the good-time intervals and within a 
region of interest (ROI) with radius of 10° from the source radio 
position applying a cut on the zenith angle parameter (< 105°) 
to avoid the Earth albedo. The following steps are to calculate the 
live-time, the exposure map and the diffuse response. 

With all these information at hands, we performed an analysis 
by using an unbinned likelihood algorithm (gt like) in three sepa- 
rate energy bands, namely 0.1-1, 1-10 and 10-100 GeV. The model 
included the isotropic and Galactic diffuse backgrounds, the source 
of interest, all the 1FGL sources in the ROI and, possibly, addi- 
tional sources not included there but identified in the map. For all 
the point sources we assumed a power law spectrum, with flux and 
photon index as a free parameter and calculated the corresponding 
test statistic (TS, see Mattox et al. 1996 for a definition; in prac- 
tice one assumes \/TS ~ a, the significance of the detection). For 
each energy bin showing a significance smaller than TS = 16 we 
derive upper limits. 

For 1ES 0229+200 and 1ES 0347-121 we only obtained up- 
per limits. For PKS 0548-322 we obtained a detection in the high 
energy (10-100 GeV) bin, upper limits in the other two bands. 

For RGB J0152+017 we have detections in all three energy 
bands (see Table [T}. In this case a fit on the whole 0.1-100 GeV 
energy band gives a spectral slope ofF = 2. 0± 0.2 consistent with 
the binned analysis described above. Since the significance of the 
detection is relatively large, we also derive a light-curve with a bin- 
size of 10 days to study possible long-term variability. The source 

2 accessible from |http : //f ermi . gsf c ■ nasa. gov| 
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is barely detected (with TS > 10) in about 10 bins. For the two 
time bins (MJD 54872-54882 and MJD 55272-55282) showing 
the highest significance, TS > 12 (~ 3.5<r), we also derived a 
spectrum, which is very different in the two cases (r = 1.15±0.2 
and F = 2.6 ± 0.4, respectively). Although the significance of 
these two spectra is rather low, their difference may suggest that 
the long-term flat spectrum is in fact the average over hard and soft 
spectra. 

In the analysis of the RGB J0152+017 data we discover sig- 
nificative (TS = 49) 7— ray emission from a source distant ~ 
1° from RGB J0152+017. The position, R.A. = 27.15°, Dec.= 
+ 1.45° (with a positional uncertainty, measured by the 95% con- 
tainment radius without systematics, of 0.06°) corresponds to that 
of the flat spectrum CRATES radio source J0148+0129, associated 
to the source PMN J0148+0129 (of unknown redshift), not included 
in the 1FGL catalogue. 

Note that the reprocessed emission component is in princi- 
ple expected to be extended (pair "echo"), while the procedure dis- 
cussed above to derive upper limits and the fluxes assumes a point- 
like source or a source with an extension contained within the LAT 
PSF. As we verified a posteriori (see i)3.2.4| ( this is a reasonable 
assumption for the cases discussed in the present work. 

2.2 Multifrequency data 

SEDs are shown in Fig.[T] We use the multifrequency data already 
considered in TlOa and the H.E.S.S. spectrum of 1ES 0548-322 
recently published by Aharonian et al. (2010). 

All the TeV spectra have been obtained by H.E.S.S. In three 
cases the sources have shown very small spectral changes dur- 
ing observations extended over several months (1ES 0347-121) or 
years (1ES 0229+200, PKS 0548-322 ). Based on this, as in TlOa 
and TlOb, we assume that the averaged H.E.S.S. spectra are repre- 
sentative of the TeV component also during the period of the LAT 
observations. 

For RGB JO 152+0 17, instead, the data have been obtained in 
a relatively short exposure (~ 2 weeks) during a ToO observation 
in 2007 (Aharonian et al. 2008), and therefore the spectrum can- 
not be considered as representative of an average state. This fact, 
together with the pronounced GeV variability suggested by LAT 
data (we report in Fig.Q]the two spectra discussed in §2.1), is prob- 
ably the reason for the apparent disagreement between the LAT 
and H.E.S.S. spectra visible in Fig.[T] Considering this caveat, for 
completeness we still use the average LAT spectrum and the sin- 
gle epoch H.E.S.S. spectrum in the following analysis although the 
inferences that we can draw about the blazar emission model and 
IGMF constraints from this source will reflect the corresponding 
uncertainties affecting the data. 

We show the most recent optical-UV and X-ray data from 
Swift, taken during or in vicinity of the first three months (August- 
October 2008) of LAT pointings. We refer the reader to TlOa for 
the details of the analysis of these data and a general discussion 
of the SEDs. For RGB J0152+017 there are no Swift observations 
close in time to the LAT ones. As in TlOa we use in the model- 
ing the X-ray spectrum close in time to the TeV observations taken 
with XMM-Newton and RXTE (Aharonian et al. 2008). 

In three cases (RGB J0152+017, 1ES 0229+200 and 1ES 
0548-322) the optical-UV data describe a rather steep contin- 
uum that, as already discussed in Tavecchio et al. (2009) for 1ES 
0229+200, traces the emission of the host galaxy. In these cases, in 
order not to overestimate the UV flux, it is necessary to assume 
a strong roll-off of the non-thermal continuum of the jet below 



the soft X-ray band. The optical to X-ray flux is very hard and 
in the most extreme case of 1ES 0229+200 requires a spectrum 
F(u) oc v 1 / 3 , as expected from the tail of the synchrotron emis- 
sion of high-energy electrons. 



3 MODELING THE SED 

3.1 Primary blazar emission model 

SEDs of TeV BL Lacs are generally adequately reproduced by 
one-zone leptonic models (e.g. Ghisellini et al. 1998, Tavecchio 
et al. 1998). Support to the idea that the emission mainly origi- 
nates in a single region comes from the existence of a characteristic 
timescales, thought to be related to the size of the emitting region^ 
and from the tight correlation between variations at very different 
frequencies (e.g. Fossati et al. 20080 

We use the one-zone leptonic emission model fully described 
in Maraschi & Tavecchio (2003). Briefly, the emitting region is a 
sphere with radius R with a tangled and uniform magnetic field B. 
Relativistic electrons are assumed to have an isotropic distribution 
and to follow a smooth broken power law energy spectrum with 
normalization K and indices n\ from 7 m in to 7b and 722 above 
the break up to 7 max . These electrons emit through synchrotron 
and synchrotron-self Compton (SSC) mechanisms. SSC emission 
is calculated assuming the full Klein-Nishina cross section. This is 
particularly important for HBLs as those considered here, in which 
the effects of the KN cross section essentially shape the SSC com- 
ponent above few tens of GeV (e.g. Tavecchio et al. 1998). Rela- 
tivistic amplification of the emission is described by the Doppler 
factor, S. These nine parameters (R, B, K, m, 712, 7mm, 7b, 7max, 
5) fully specify the model. 

Once the intrinsic spectrum F[ n t(E) is derived, the cor- 
responding observed spectrum is calculated as Fi nt , bs(E) = 
F int (E)e- T(E) , where r(E) is the energy-dependent optical depth 
(specified below) for the 77 — > e interaction with the optical-IR 
photons of the extragalactic background light (EBL). 

3.2 The reprocessed spectrum 

3.2.1 Assumptions 

During their travel to the Earth, a fraction of the blazar primary 
high-energy 7-rays interacts with the low energy (optical-IR) ex- 
tragalactic photons producing electron-positron pairs. The pairs, in 
turn, inverse-Compton scatter the photons of the CMB, emitting 
a reprocessed high energy component. The flux and the spectral 
shape of this component depend on the spectrum of the primary 
7-rays and on the intensity and geometry of the intergalactic mag- 
netic field which, deviating the pairs during their cooling, affects in 
an energy-dependent way the observed flux. The comparison of the 
GeV data with the expected reprocessed spectrum provides a way 
to constrain the IGMF. We take into account these effects refining 
the simplified treatment of T 10b. 

The interaction of a primary photon with energy (measured 
in TeV) EtcV with a low frequency (optical-IR) photon of the 
(EBL) results in an electron positron pair with Lorentz factor 
7 ~ E/2m e c 2 ~ 10 6 -E>rcV. These pairs, in turn, will inverse 

3 but there are exceptions, e.g. Mkn 501 (Albert et al. 2007) and PKS 2155- 
304 (Aharonian et al. 2007) showing event with very short timescales. 

4 exceptions are the so-called orphan flares (e.g. Krawczynski et al. 2004). 
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Figure 1. Spectral Energy Distributions of the sources discussed in this work. Open squares are historical data (green in the electronic version), filled triangles 
reports Swf/f/UVOT and XRT data and filled circles (red) shows the TeV spectra. The open (green) circles in the SED of RGB J0152+017 report the X-ray 
spectrum taken in the vicinity of the H.E.S.S. detection. Fermi/LAT data are those derived by our analysis. For RGB J0152+017 we also show (gray bow- 
ties) two spectra accumulated over ten days (54872-54882, dotted; 55272-55282, dashed-dotted) showing huge variations both in flux and slope (see text). 
See Tavecchio et al. (2010a) for references. The TeV spectrum of 1ES 0548-322 is from Aharonian et al. (2010). The solid line shows the emission model 
reproducing the observed SED, including synchrotron, SSC and reprocessed emission. The dashed line reports the intrinsic SSC component. See text for 
details. 



Compton scatter the photons of the CMB producing 7-rays of en- 
ergy e ~ 7 2 /ii/ CM B ^ 2.8 /cTcmb 7 2 = 0.63 £| cV GeV. 

In TlOb we described the reprocessed 7-ray spectrum as a 
power-law (with slope 0.5 in energy) with an abrupt cut-off above 
an energy e ma x- This approximation was motivated by the fact that 
this component derives from electrons which have a very short 
cooling timescale and therefore quickly form a completely cooled 
distribution. Although this approximation was appropriate for the 



purposes of TlOb, here we adopt a more self-consistent approach, 
in order to have a better description of the reprocessed spectrum 
around e max . 

Another important aspect to consider is the possibility to have 
more than one generation of pairs. In TlOb we assumed that the 
maximum energy of the first generation reprocessed photons, e max , 
was below the threshold for further interaction with the EBL, which 
for sources at z ~ 0.1 corresponds to e max ~ 300 — 400 GeV. This 
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translates into a maximum energy of the intrinsic blazar emission, 
i5 max ~ 20 TeV. When realistic SSC spectra are considered this 
condition could be violated, since there could be an important tail 
of the intrinsic SSC component above this energy. Therefore it is 
necessary to consider that also the 7— ray photons produced through 
reprocessing can be further absorbed by the interaction with the 
EBL, producing a second generation of inverse Compton emitting 
pairs. Again, this second component could be absorbed and so on, 
possibly leading to the development an electromagnetic cascade. 

A complete, self-consistent treatment of the cascades could be 
done with Monte Carlo simulations (e.g. Coppi & Aharonian 1997, 
Elyiv et al. 2009, Dolag et al. 2009, 2011, Neronov et al. 2010, 
Taylor et al. 2011). However for sources at the redshift around 
0.1, as those considered here, a real electromagnetic cascade (i.e. 
more than two generations) would develop only if the maximum 
energy of the primary SSC radiation is very large, i? m ax > 100 
TeV. This can be derived from the relation between the energy 
of the primary, EtcV (measured in TeV) and of the reprocessed 
e photon e ~ 0.63i5xeV GeV. We make the simplified assump- 
tion that the absorption of the radiation produced by pairs becomes 
important only if its maximum energy is above a threshold eth 
at which the optical depth is larger than one. Therefore an im- 
portant second generation of pairs will develop only if the maxi- 
mum energy of the photons produced by the first-generation satis- 
fies £ m ax = 0.63-Emax > e th, namely, if the maximum energy of 
the primary photons exceeds: £ max > (e t h/0.63) 1/2 (again, E is 
in TeV, e in GeV). Generalizing this expression, the condition for 
the effective absorption of the radiation of the n-th generation of 
pairs (and therefore for the production of an important population 
of pairs of the generation n + 1) is: 



-En 



> 1600 



eth 



1.6 x 10 6 



TeV 



(1) 



Assuming e t h = 500 — 1000 GeV, a value suitable to describe the 
case of the sources analized here (see below), to have more than two 
generations, n — 2 in Eq.[T] -E max must be larger than ~ 250 TeV 
(due to the small exponent this value is not strongly dependent on 
e t h, the assumed threshold for important absorption). In the context 
of the standard SSC model used below such large energies are very 
difficult to achieve: therefore we can limit the calculation of the 
reprocessed spectrum to the second generation. 

A further assumption concerns the variability of the primary 
TeV emission. As noted by Dermer et al. (2010), the reprocessed 
radiation in the GeV range reaches the observer with delays with 
respect to the direct unabsorbed radiation. With magnetic fields of 
the order of 10 -15 G such as those derived below, the delays are 
of the order of some million years. Therefore, to compare the flux 
of the reprocessed component with the observer primary flux we 
have to assume that the primary 7-ray flux is on average stable 
over times of the order of 10 7 years. Smaller times would result in 
correspondingly lower values for the IGMF. For the case of 1ES 
0229+200, Dermer et al. (2010), considering a lower limit of t > 1 
year for the activity timescale of the source findB > 3 x 10~ 19 G. 

In the calculation the use of a specific model for the EBL is 
unavoidable. In the range of wavelength interesting here, 0.1 — 15 
fim, several of recent EBL estimates converge: the recent models of 
Franceschini et al. (2008), Gilmore et al. (2009), Finke et al. (2010), 
Dominguez et al. (2010), agree well (but see Orr et al. 2011) and 
are consistent with the low level of the EBL suggested by the recent 
observations of Cherenkov telescopes (e.g. Aharonian et al. 2006, 
Mazin & Raue 2007; see also Kneiske & Dole 2010). We adopt the 
LowSFR model of Kneiske et al. (2004) which provides an optical 



depth similar to that of all the other updated models up to energies 
of 4-5 TeV. For larger energies the predicted optical depth is lower 
than that of the other models. In this respect in our calculation the 
importance of the reprocessed component is minimized and thus 
our results can be considered as conservative. 

Finally, we note that in our derivation we are implicitly as- 
suming that the magnetic field is oriented perpendicularly to the 
direction of the relativistic pairs and maintains its coherence within 
the region in which pairs emit and cool. A non perpendicular mag- 
netic field is less effective in deviating the emitting pairs, therefore 
determining a higher final observed flux. On the other hand, if the 
coherence length As of the field (namely the size of the regions 
in which the field maintain its coherence) is smaller than the typi- 
cal cooling length of the pairs ct coo i ~ O^^ 1 Mpc, the resulting 
magnetic field will be correspondingly larger, being in this case 
B oc \ B (e.g. Neronov & Semikoz 2009). In any case, both ef- 
fects would result in a larger value for the IGMF intensity making 
our lower limits conservative estimates. 



3.2.2 First generation spectrum 

We calculate the energy distribution of the emitting pairs assum- 
ing the equilibrium (ensured by the very small cooling time of the 
pairs) between injection of new energetic pairs in the volume be- 
tween the source and the Earth (produced through the conversion 
of absorbed 7-rays) and their cooling through IC emission. From 
the general kinetic equation for the electron distribution N(^/) we 
get: 



iV( 7 ) = 



7 



(2) 



where the injection term Q(7) is calculated considering the absorp- 
tion process: 



Qd) 



E 



-t(E) 



with E = 2m e c 7 (3) 



where Fi n t(E) is the spectrum of the intrinsic (primary) blazar 
emission and t(E) is the optical depth for absorption of 7-rays. 
The cooling term 7 in Eq.[2]is specified assuming that the cool- 
ing is dominated by the IC scattering on the CMB photons in the 
Thomson limit: 



7 ! 



4 <tt tt 2 

UCMB7 • 

3 m c c 



(4) 



where C/cmb is the CMB energy density at the redshift of the re- 
processing region. 

Once the pair energy distribution is specified through Eq. [2] 
we calculate the corresponding inverse Compton spectrum Fi(e). 
The normalization of the spectrum is derived considering that, as 
said above, the system is in the so-called fast cooling regime, i.e. 
the total (i.e. energy integrated) emitted reprocessed power is equal 
to the total power injected into pairs (neglecting the rest mass of 
pairs), in turn equal to the energy integrated power of the absorbed 
radiation (e.g. TlOb). 

Note that since the primary blazar spectrum Fi n t(E) is usu- 
ally hard in the range of energies where absorption is important 
and t(E) increases rapidly with energy, the spectrum of injected 
electrons will also be hard and therefore, consistently with the as- 
sumption of TlOb, the electron energy distribution resulting from 
Eq. (O will be a simple power law ^(7) oc 7~ 2 (since the integral 
is almost independent on the lower limit 7), except for the high- 
est energies, where the distribution will have a rapidly decreasing 
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7min 7b 7max ni "2 B [G] K [cm 3 ] _R [cm] 6 -Bigmf [G] 
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7.3 x 10 8 


5.8 x 10 16 


34 


7 


x 10" 


14 


1ES 0347-121 


5 x 10 4 


1.5 x 10 6 


2 x 10 6 


2.2 


3.1 


2 x 10" 2 


1.8 x 10 4 


3.35 x 10 16 


15 




10 -14 


PKS 0548-322 


1.5 x 10 4 


2.8 x 10 5 


1.5 x 10 6 


2.0 


4.1 


0.1 


10 4 


3.6 x 10 15 


16 


5 


x 10" 


15 



Table 2. Input parameters for the emission models shown in Fig[T] See text for definitions. For 1ES 0229+200 we assume a simple power law electron energy 
distribution (see Tavecchio et al. 2009) and we also report (B) the parameters used for the model reproduced in Fig.[6]and discussed in the text. 
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Figure 2. Ratio of the 7-ray mean free path for photons of E = 10 TeV 
and source distance versus redshift using the LowSFR model of Kneiske et 
al. (2004). Stars indicate the redshift of the sources considered in our work. 



tail. Importantly, this also implies that the number of pairs from 
the conversion of the primary photons at the highest energy, i? max , 
exceeds the number of pairs deriving from photons af lower en- 
ergies. Therefore, the reprocessed spectrum is dominated by the 
emission from the pairs originated by the absorption of the primary 
photons at i5 m ax, allowing to identify a characteristic distance from 
the source (<f 7 , the mean free path for photons of energy i5 max ), at 
which the entire first generation reprocessed spectrum is produced. 

The effect of the IGMF is to curve the trajectory of the emit- 
ting pairs, resulting in a dilution of the reprocessed spectrum within 
a solid angle which is a function of the pair energy and IGMF inten- 
sity, Sl 7 = 27r(l — cos 8-y). The angle (9 7 can be estimated assuming 
it is the angle by which the electron velocity vector changes in the 
cooling length rf coo i: 6» 7 ~ ct cool /r L = 1.17 B-is(l + z r )~ A ^ 2 
where vl = jrrieC 2 / (eB) ~ 2 x 10 24 76-BZi 5 cm is the Larmor 
radius of the electron. 

In TlOb we assume that for values of the energy and IGMF 
for which the angle of curvature due to the IGMF is within the 
beaming solid angle of the intrinsic blazar emission, Q c ~ ir8 2 , the 
observed spectrum is not affected by IGMF. Conversely, if the pairs 
are completely isotropized during their lifetime the observed flux 
will be suppressed by a factor n c /47r. A convenient approximation 



for the observed first generation reprocessed spectrum, including 
both the absorption by the EBL and the effect of the IGMF on the 
emitting pairs, is calculated from F\ (e) as: 



F lt0hs {e) = Fi(e) 



-ri(e) 



n c + fi 7 (e) 



(5) 



where, Fi(e) is the IC spectrum calculated above and, again, Q c 
is the solid angle of the primary blazar emission and f2 7 (e) is the 
solid angle (dependent on the photon energy and IGMF intensity) 
into which the radiated energy of a pair is spread (see TlOb for 
details). For the models reported here we assume 8 C — 0. 1. Smaller 
values of the beaming angle would result in lower fluxes for the 
fully isotropised reprocessed component. 

As long as the pairs of the first generation are produced and 
cool at distances from the blazar small compared to the distance 
of the blazar to Earth (d 7 /D ~ 1/t <§C 1, where d 1 is the mean 
free path of photons as measured from the blazar and D the blazar 
distance) we can make the approximation ri ~ r. Of course the 
validity of this approximation depends on the considered energy 
and it will be violated at low energy, for which the optical depth for 
the primary photons is t(E) < 1. However, as discussed above, we 
can assume that the reprocessed emission mainly derives from pairs 
produced by the conversion of primary 7-rays at i5 m ax, which is 
typically around 10 TeV. For illustration, in Fig.[2]we report the ra- 
tio d-f/D for photons of energy E = 10 TeV for the adopted EBL 
model for the range of redshift interesting for the sources consid- 
ered here. For 1ES 0229+200 and 1ES 0347-121, (z = 0.14 and 
0.18, respectively) d-,/D ~ 0.1 — 0.2, ensuring that the approxima- 
tion t\ m t is good. For the other two sources (z = 0.07 and 0.08), 
instead, the approximation is not applicable. However, for the same 
reason - the optical depth is low - only a small fraction of the pri- 
mary continuum is reprocessed and therefore the contribution of 
the reprocessed emission to the total observed GeV-TeV spectrum 
is not important (see Fig. [3] the maximum total contribution of the 
reprocessed component is less than 30%). 

We finally note that the approximation considering the blazar 
emission pattern as a uniform cone is rather useful in simplifying 
the calculation and provides correct results when we are consider- 
ing integrated quantities. However, the right pattern should be con- 
sidered when one is interested in the detailed spatial distribution of 
the reprocessed component (see ^3.2.4t . 



3.2.3 Second generation spectrum 

The observed spectrum emitted by the second generation of pairs, 
i 7 2,obs(e), is calculated iterating the procedure described above for 
the first generation. 

A further assumption that we use is that the emitting pairs of 
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Figure 3. High energy SED of the sources discussed in this work. Open triangles are the observed TeV data, open squares are the data-points corrected 
for absorption by the EBL assuming a distance equal to that of the source. Diamonds are the LAT data-points derived by our analysis: the horizontal lines 
indicate the band considered to derive the flux. The solid thick line (black in the electronic version, labelled "Total") shows the emission model reproducing 
the observed SED, including SSC and reprocessed emission. Thin orange lines shows the intrinsic (dashed, A) and observed (solid, B) SSC spectrum, the 
thin blue lines are the intrinsic (dashed, C) and observed (solid, D) first generation reprocessed spectrum. The dotted line is the first generation reprocessed 
spectrum assuming no IGMF. The green thin line (E) is the second generation reprocessed spectrum, providing negligible contribution in all the cases. See text 
for details. 



the second generation are initially (i.e. when created) still colli- 
mated within the beam of the primary blazar emission. This re- 
quirement is generally satisfied since the absorbed first-generation 
7-rays converting into the second generation of emitting pairs are 
highly energetic, e > 100 — 200 GeV, and they are produced by en- 
ergetic (7 > 10 r ), fastly-cooling pairs of the first generation. These 
pairs, during their short lifetime are only slightly deviated from the 
original direction which was that of the primary blazar emission, 
the deviation angle being 9 < 0.6B-isj^ 2 deg. Only for values 
of the IGMF exceeding those considered in the following, the tra- 
jectories of the first generation pairs at the highest energy will be 
substantially curved during their lifetime and this approximation 
would break down. 



3.2.4 Angular extension 

The spreading due to IGMF implies that the reprocessed emission 
is not point-like but it is characterized by a finite angular size (e.g. 
Elyiv et al. 2009) that, especially at the higher energies (above 10 



GeV), where the angular resolution of LAT is better (with contain- 
ment radius of few tenths of degree), can exceed the PSF width. Un- 
der these conditions the reprocessed component should be treated 
as an extended source and the standard LAT analysis used above to 
derive upper limits or fluxes cannot be applied. 

It is therefore important to quantify the expected angular 
broadening as a function of photon energy and IGMF intensity and 
compare it with the LAT PSP0. Under the simplifying assumption 
that the intensity of the primary emission is uniform within the cone 
with semi-aperture 6 C (see Hig. [4j a simple relation for the typical 
angular size of the reprocessed component, 9 V can be derived (see 
also Neronov & Semikoz 2009): 



9v(e) 



^7 

1' r 



(6) 



where we expressed the ratio between the source distance and the 
typical interaction length of photons with energy E max (see jmj 



|http: //www-glast . slac . stanf ord. edu/software/IS/glast_Lat-perf on 
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as t ~ D/d 7 > 1 (see Fig.[4]l. Eq.[6]simply states that for energies 
and IGMF intensity for which the pair deviation is larger than the 
cone semi-aperture we can see the entire conversion surface, that is 
the spherical shell described by the angle 9 C projected into the sky. 
On the contrary, when the deflection is less than 9 C , we only see the 
(energy dependent) portion of the surface depending by the angle 

This simple estimate is derived under the assumption that the 
blazar uniformly illuminates the associated beaming cone. How- 
ever this is a crude approximation: indeed the blazar relativistically 
boosted intrinsic emission is strongly anisotropic, being concen- 
trated along the jet axis. To properly estimate the expected bright- 
ness profile of the reprocessed emission one should then consider in 
detail the emission pattern, that for the integrated flux is described 
by <5 4 , where S is the beaming Doppler factor. Due to the strong an- 
gular dependence of 5 the emission pattern of the intrinsic emission 
is rather "narrow" and the observed extension size will be in gen- 
eral less than what derived through Eq|6] In the following we will 
derive the expected profile under the simplifying assumption that 
the observer is exactly on-axis. More general cases will be consid- 
ered elsewhere. 

The flux of the reprocessed component as measured by the 
observer from an annulus between 9 and 9 + d9 can be written as: 

f{e,8)d8 = 2irsin8I(e,8)d8 (7) 

for angles 9 < d 1 (e)/r (for larger angles the observer does not 
receive flux since the pairs are not deflected enough to point to- 
ward the observer). The problem reduces then to calculate 9), 
the spectrum of the reprocessed component at a given observation 
angle. Analogously to the procedure discussed above (for an ob- 
served on axis), 7(e, 9), is derived from the corresponding primary 
intensity at the angle 9' , I(E, 9'). The relation between the two an- 
gles, 8' (as measured from the source) and 9 (as measured form the 
observer) can be derived by geometrical arguments, 8' ~ (r — 1)8 
assuming, as discussed above, that the reprocessed emission orig- 
inates at a typical distance d 7 . In turn, I(E, 9 1 ) is calculated from 
that measured on-axis, I(E, 0) using the standard beaming rela- 
tions. 

Having derived /(e, 8), to evaluate the spatial size of the ex- 
tended halo we calculated the (energy dependent) angle within 
which one collects half of the entire flux, 80.5(e) according to the 
following condition: 

r8 05 (c) roc 

/ f(e,0)d0 = 0.5/ f{e,8)d8 (8) 

Jo Jo 

Curves reporting 80. 5(e) for different values of T and Bigmf 
for 1ES 0229+200 and 1ES 0347-121 are compared with the LAT 
containment radius at 68% and 95% in Fig. [5] (we do not con- 
sider the other two sources since, in our models, the GeV flux is 
dominated by the point-like blazar emission). Not surprisingly, the 
curves display the same behavior of the simple estimate given by 
Eq.[6] At low energy #0.5 does not depend on the energy because 
the emission of the corresponding electrons is isotropic and the ex- 
tension of the reprocessed emission is dictated by the typical beam- 
ing angle of the primary continuum, ~ 1 /T. At sufficiently large 
energy, instead, the electrons start to emit within small angles and 
thus the extension is fixed by the corresponding deviation angle 
of the electrons, # 7 oc e _1 . The values of #0.5 are always within 
the containment radius of LAT and therefore we can safely use the 
standard analysis to derive the upper limits. 



4 RESULTS 

The models used to reproduce the observed SEDs are reported in 
Fig. E and the input parameters are listed in Table 1 . As detailed 
in §3.2, in reproducing the high-energy part of the SEDs we con- 
sider the contribution of the absorbed intrinsic SSC component, 
that of the absorbed first generation reprocessed component and we 
also consider the second-generation emission (although providing 
a negligible contribution in most cases). The latter two contribu- 
tions depend on the assumed SSC spectrum and on the value of the 
IGMF intensity, Bigmf (also reported in Tab. 1). For simplicity, in 
the SEDs we only report the intrinsic SSC component (dashed line) 
and the total (intrinsic plus reprocessed) observed spectrum (solid 
line). A zoom on the high-energy band, reporting all the different 
components, is reported in Fig.[3] showing the intrinsic SSC spec- 
trum (orange), the first generation reprocessed spectrum (blue) and 
the second generation spectrum (red). In all the cases the dashed 
lines report the emitted spectrum, the solid line the spectrum af- 
ter absorption by interaction with the EBL. The dotted line cor- 
respond to the first-generation reprocessed spectrum in absence of 
any IGMF. The black thick solid line is the total observed spectrum. 
For comparison in Fig.[3]we report both the observed TeV data and 
the data deabsorbed assuming the Kneiske et al. (2004) LowSFR 
model. 

Even if possibly overestimated with our assumptions, the con- 
tribution of the second-generation component is almost negligible 
in all the cases (but for 1ES 0229+200). Conversely, the first gen- 
eration component is generally important, especially at energies 
around few hundreds GeV. The importance of the contribution of 
this component mainly depends on the shape and maximum energy 
of the intrinsic SSC emission: a hard spectrum peaking at large 
energies (10 TeV) will result in a large absorbed power and, conse- 
quently, a large power in the reprocessed component which extends 
in the TeV band. On the other hand, if the primary SSC spectrum 
peaks below 10 TeV the reprocessed component will contribute 
mainly at sub TeV energies. Generally for values of the IGMF 
above Bigmf = 10~ 15 G the contribution of the reprocessed emis- 
sion in the 0.1-100 GeV band goes below the present sensitivity of 
LAT. 

The most stringent lower limit to the IGMF intensity is pro- 
vided by 1ES 0347-121 (that located at the largest distance): in or- 
der not to exceed the upper limit in the 10-100 GeV band one has to 
use an IGMF intensity of Bigmf = 10~ 14 G. For 1ES 0229+200 
the same energy bin constrains Bigmf = 2 x 10" 15 G (compara- 
ble to that derived in TlOb). For the other two sources, due to the 
smaller distance the contribution of the reprocessed component is 
less important and basically the value of B is unconstrained. 

Since the level of the reprocessed component depends on the 
primary SSC spectrum, also the limits on the IGMF based on the 
contribution of this component to the overall spectrum will depend 
on it. As an example of the importance of the intrinsic SSC spec- 
trum in shaping the reprocessed component and of the uncertain- 
ties in the derived parameters we report in Fig. [6] an alternative 
fit for 1ES 0229+200 assuming a SSC component peaking at 10 
TeV. The first generation reprocessed component extends well in 
the TeV band, dominating the spectrum below 5 TeV. Also the sec- 
ond generation component is now very powerful with a flux of the 
same order of that of the first generation component. In order not to 
exceed the level of the data around 200 GeV one is forced to assume 
a very large IGMF, B — 7 x 10" 13 G. We stress that in this case 
the limit to the IGMF intensity is constrained by the data-points at 
100-200 GeV, not by the LAT upper limits. 
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Figure 4. Sketch of the geometry assumed for the calculation of the spatial extension of the reprocessed emission. See text for definitions. 




Figure 5. Expected size of the reprocessed emission (expressed as the angle containing half of the total flux, 9 0.5) as a function of photon energy for 
1ES 0229+200 and 1ES0347-121 for different values of the blazar bulk Lorentz factor T and IGMF intensity. For comparison also the 68% and 95& LAT 
containment radius are reported. 



This example makes clear how the derivation of the intrinsic 
spectrum of the source (and hence of the physical parameters of 
the emitting region) and the estimate of the IGMF are interlaced. In 
particular it makes clear that some of the models presented in TlOa 
(in particular for the cases of 1ES 0229+200 and 1ES 0347-121) 
are no longer adequate to reproduce the data once the reprocessed 
component is taken into account, unless one assumes large values 
for the IGMF. As in TlOb, here we prefer to adopt a conservative 
approach, trying to minimize the value of the required value of the 
IGMF. Therefore, the model reported in Figs. [TJ and [5] assumes a 
SSC peak at energies lower than assumed in TlOa. This choice, 
in turn, implies that we derive value of the magnetic field in the 
emission region larger than that found in TlOa. 

In Fig. [7] we report the updated plot from TlOa showing the 



value of the magnetic field in the emission region and of the break 
Lorentz factor of the relativistic electrons (those emitting at the 
SED peaks) for the sample of the high-energy emitting BL Lacs 
considered in TlOa. For the four sources considered in the present 
work we show both the values derived in TlOa and those derived 
here. The different symbols indicate: sources detected both by LAT 
(belonging to the LBAS list Abdo et al. 2009) and Cherenkov tele- 
scopes (red open circles), by LAT only (open green squares) and 
by Cherenkov telescopes only (i.e. the four sources considered in 
this work, triangles). For the latter sources we report both the val- 
ues derived in TlOa (filled triangles) and those inferred here (open 
triangles). As discussed in TlOa, the sources of the latter group 
populate an extreme region of the plane B — 7b, at very low values 
of the magnetic field, B < 10~ 2 G, and high electron energies, 
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Figure 6. High-energy SEDs of 1ES 0229+200 reproduced with a model as- 
suming a large peak energy (10 TeV) for the primary SSC component. The 
consequent high level of the reprocessed emission (dominating the overall 
emission at TeV energies) requires a large value of the IGMF to not exceed 
the observational points below 5 TeV. 



7b > 1CP J . As discussed above, in order to minimize the impor- 
tance of reprocessed component we assume for 1ES 0229+200 and 
1ES 0347-121a SSC component peaking at lower energies than 
those adopted in TlOa. This, in turn means a larger magnetic field, 
as shown by the position of the new points in the plot. For RGB 
0152+017 and PKS 0548-322, instead, the small redshift makes the 
reprocessed component less important than in the other two sources 
and we can model their SEDs with the same parameters of T 10a. 
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Figure 7. Magnetic field versus the break Lorentz factor 7b for the high- 
energy detected BL Lacs (LBAS and Cherenkov) modeled with the one- 
zone synchrotron-SSC model, updated from Tavecchio et al (2010b). Red 
open circles (triangles) show the values for the known TeV sources detected 
(not detected) by LAT. Open green squares are for the LAT BL Lacs not de- 
tected in the TeV band. For 1ES 0229+200 and 1ES 0347-121 we report 
both the points corresponding to the parameters assumed in T 10a (filled tri- 
angles) and those derived here (open triangles) considering the reprocessed 
component. For the other, close-by sources, the parameters adopted in TlOa 
adequately reproduce the SED even with reprocessing. 



5 CONCLUSIONS 

The inclusion of the emission arising from the reprocessing of the 
absorbed primary 7 rays can play an important role in the detailed 
modeling of the SEDs of blazars. In particular we have shown that 
for sources with hard primary VHE continuum and redshift around 
0. 1 the reprocessed component can even dominate below few TeV, 
unless the value of the IGMF intensity exceeds 10~ G. 

The possibility that the sub-TeV radiation from blazars re- 
ceives an important contribution for the reprocessing of the high- 
energy power absorbed through interaction of high-energy photons 
with the infrared-optical background has a number of important 
consequences. The first obvious fact is that the reprocessed com- 
ponent should be considered in the modeling. However, since the 
exact shape of this component is linked to the unknown value of 
the IGMF and the often not well determined maximum energy of 
the primary SSC component, this introduces several unknown pa- 
rameters in the fit, breaking the simplicity of the one-zone leptonic 
model and preventing to obtain a unique set of physical parameters. 

We would like to stress that to constrain the IGMF using the 
limits to the reprocessed component it is mandatory to consider 
realistic SSC spectra for the primary radiation. In particular, the 
widely considered one-zone leptonic model with typical parame- 
ters hardly predicts SSC peaks above 20-30 TeV, even in the case 
of the most extreme HBL objects as those considered here. This 
fact limits the total power that can be reprocessed into the GeV 
band and implies that, for the sources located at typical distance 
z — 0.1 — 0.2, a real electromagnetic cascade does not develop 
and the reprocessing process ends after the second generation of 



pairs. However, blazar primary spectra could extend well above the 
limits derived in the SSC model discussed here if, for instance, the 
VHE emission is the result of hadronic processes (e.g. Miicke et 
al. 2003). In this case, the reprocessing could extend beyond the 
second generation, making the use of Monte Carlo calculations un- 
avoidable. In this case, if the primary spectrum extends well above 
10 TeV with a hard slope, the total amount of reprocessed radiation 
increases, implying larger values of Bigmf- 

The fact that the reprocessed component can provide an im- 
portant contribution even at energies in the TeV range can have 
a number of important consequences for the next generation of 
planned Cherenkov facilities (CTA0). In particular, this fact high- 
lights the possibility to obtain improved limits to the intensity of 
the IGMF by using precise measures of the overall VHE spectrum 
alone, especially if low thresholds can be reached. 

A natural and testable consequence of the idea that the emis- 
sion of blazar above 50-100 GeV is a mixture of the primary (at 
the highest energy) and reprocessed (at lower energies) emission is 
the possibility to have spectra more complex than the widely as- 
sumed simple power law, with even multiple bumps marking the 
peak of the emission of different generations of pairs. Variability 
could also be used to disentangle different components in the 7-ray 
spectrum: indeed, as long as the magnetic field exceeds B = 10 -20 
G, the emitting pairs are effectively deviated and, consequently, the 
reprocessed emission is diluted over times that can reach w 10 7 
years around 1-10 GeV (e.g. Neronov & Semikoz 2009, Dermer 
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et al. 2010). Therefore we expect that the reprocessed emission is 
rather stable and can be easily distinguished from the highly vari- 
able blazar primary continuum. 
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